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ABSTRACT 

X-ray emission line profiles provide the most direct insight into the dynamics 
and spatial distribution of the hot, X-ray-emitting plasma above the surfaces of 
OB stars. The O supergiant ( Puppis shows broad, blueshifted, and asymmetric 
line profiles, generally consistent with the wind-shock picture of OB star X-ray 
production. We model the profiles of eight lines in the Chandra HETGS spec¬ 
trum of this prototypical hot star. The fitted lines indicate that the plasma is 
distributed throughout the wind starting close to the photosphere, that there is 
significantly less attenuation of the X-rays by the overlying wind than is gen¬ 
erally supposed, and that there is not a strong trend in wind absorption with 
wavelength. 

Subject headings: line: profiles — stars: early-type — stars: mass loss — stars: 
winds, outflow — stars: individual (£ Puppis) — X-rays: stars 

1. Introduction 

The nature of the copious soft X-ray emission from hot stars has been a longstanding 
controversy since its discovery in the late 1970s (Harnden et al. 1979; Cassinelli & Olson 
1979). Solar-type coronal emission was first assumed (Cassinelli & Olson 1979; Waldron 
1984), but wind-shock models of various types gained currency throughout the following 
decade (Waldron 1984; Owocki, Castor, & Rybicki 1988; MacFarlane & Cassinelli 1989; 
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Chen & White 1991; Hillier et al. 1993; Cohen et al. 1996; Feldmeier et al. 1997; Feldmeier, 
Puls, & Pauldrach 1997). More recently, hybrid magnetic wind models have been proposed 
for some hot stars (Gagne et ah 1997; Babel & Montmeric 1997; ud-Doula & Owocki 2002). 

Until the launch of Chandra and XMM, no observational diagnostics were available that 
could provide a direct discriminant between the coronal and wind-shock paradigms. But due 
to their superior spectral resolution, both of these new telescopes allow for the separation of 
individual emission lines and the resolution of Doppler-broadened line profiles. The spectral 
resolution of Chandra’s grating spectrometers exceeds -A- ~ 1000 (for the FWHM) which 
corresponds to a velocity of 300 km s -1 , and that of the XMM RGS is almost as great. This 
compares favorably to the terminal velocities of the radiation-driven winds of O stars, which 
approach v <*, = 3000 km s -1 , implying ~ 20 resolution elements for a velocity range of 2^00. 

At the most basic level, the X-ray emission lines from hot stars will either be narrow 
and therefore roughly consistent with coronal emission or broad and roughly consistent with 
wind-shock emission. This is because in the wind-shock model the high velocities of X-ray- 
emitting plasma embedded in the wind would Doppler-shift the emission across a range of 
wavelengths. Initial papers reporting on Chandra and XMM observations of various O and 
B stars (Schulz et al. 2001; Kahn et al. 2001; Waldron & Cassinclli 2001; Cassinclli et al. 
2001; Miller et al. 2002; Cohen et al. 2003) discussed line widths, which vary from large for 
early O stars to small for B stars. Some of these initial studies noted that the emission 
lines can be blueshifted and the profiles somewhat asymmetric, but none of these studies 
discussed or modeled the shapes of the resolved emission lines. 

In this paper we fit a specific model of X-ray emission line profiles in an expanding, 
emitting, and absorbing wind to a Chandra HETGS/MEG spectrum of Q Puppis. The 
model we fit is empirical and flexible, with only three free parameters. The model assumes 
a two-component fluid, having as its major constituent the cold, X-ray-absorbing plasma 
that gives rise to the characteristic UV absorption lines observed in hot star winds, and 
as its minor constituent the hot, X-ray-emitting plasma. This empirical model is not tied 
to any one specific physical model of X-ray production, and is general enough to fit data 
representative of any of the major models, so long as they are spherically symmetric. To 
the extent that wind-shock models are found to be consistent with the observed line profiles, 
our model parameters can be used to constrain the physical properties of the shock-heated 
plasma. This ultimately can be used to constrain the values of the physical parameters of 
the appropriate wind-shock model. 

In section 2 we discuss the physical effects leading to non-trivial line shapes and describe 
the specific empirical model we use to perform the fits. In section 3 we describe the Chandra 
( Puppis dataset and how we perform the line fitting and parameter estimation. And in 
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section 4 we discuss the derived model parameters in the context of the various physical 
models that have been proposed to explain hot-star X-ray emission, as well as in the context 
of other X-ray diagnostics that have recently been applied to the data from this star. 


2. Theoretical Considerations 

In the context of the fast, radiation-driven winds of OB stars, a source of X-ray emission 
embedded in the wind will lead to Doppler broadened profiles, but only the X-ray emitting 
plasma traveling at the wind terminal velocity directly toward or away from the observer 
will lead to maximal blueshifts and redshifts. The amount of emission at each intermediate 
wavelength, and thus the shape and characteristic width of the line, depends on the spatial 
and velocity distribution of the hot plasma. As described by MacFarlane et al. (1991) in 
the case of a shell of X-ray emitting plasma in a hot star wind, the continuum absorption 
of X-rays by the cool component of the wind will cause the resulting emission lines to be 
attenuated on the red side and be relatively unaffected on the blue side of the line profile. 
The apparent peak of the emission line thus shifts to the blue side of line center, and the 
line is asymmetric, with a shallower red wing and a steeper blue wing. 

This basic idea of Doppler broadened emission from a hot wind component and contin¬ 
uum absorption by the cool wind component leading to a broadened, shifted, and asymmetric 
line was extended from a shell to a spherically symmetric wind by Ignace (2001). He showed 
that model line profiles could be generated analytically for a constant-velocity wind. Owocki 
& Cohen (2001) extended this concept further, to an accelerating wind, with a model having 
four free parameters. Two describe the spatial distribution of the X-ray emitting plasma - 
R 0 , the minimum radius of X-ray emission and q, the radial power-law index of the emissiv- 
ity. There is assumed to be no emission below r = R 0 . Above r = R 0l volume emissivity 
is assumed to scale like the density of the wind squared p 2 (since collisional processes and 
recombination dominate the ionization/excitation kinematics), with an extra factor r~ q al¬ 
lowing for spatial variation of shock temperatures, cooling structures, and density and filling 
factor of the shocked material. The parameter (3 controls the velocity of the wind, which is 
assumed to follow a “beta-velocity law”: 

v(r) = v^l - Rjr)' 3 . (1) 

Both components of the wind follow the same velocity law in all cases discussed here, but in 
principle they could be allowed to differ. The fourth parameter, r*, characterizes the amount 
of absorption in the wind 
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where k is the line opacity or mass absorption coefficient (cm 2 s x ). This relates to the 
commonly-quoted radius of optical depth unity Ri by the equation (for (3 = 1) 


R 1 

-R* 


1 

1 - exp (-1/r*) 


~ r* + 0.5 for r* > 0.5. 


( 3 ) 


The error in the above approximation is less than 10% for r* > 0.67 (or R\/R* > 1.3). 

The model put forward by Owocki & Cohen (2001) is a phenomenological one. It 
describes the physical properties of the hot and cool components of the wind, but does 
not describe the physics underlying the generation of the hot plasma. It is therefore quite 
flexible, capable of describing a thin shell of X-ray emitting plasma, including a coronal- 
type zone near the photosphere, as well as wind shocks distributed spatially throughout the 
wind, with the shock distribution and wind velocity varying with radius. This model is 
thus capable of constraining properties of both the X-ray-emitting and the X-ray-absorbing 
wind components. When applied to an ensemble of lines it has the potential to constrain 
these properties as a function of both temperature and wavelength. Assuming turbulent 
and thermal broadening are negligible and combining the other model assumptions, the line 
profile as a function of scaled wavelength x = (A/Ao — l)(c/uoo) is given by 

r °° r —(9+2) 

Lx(X 7— exp [ -r (/^> r )] dr > ( 4 ) 

J r=r x %J 

where r x = max [R 0 ,R*/( 1 — |x| p// ^)], fi x = x/(l — _R*/r) /3 and r(fi,r) (which is propor¬ 
tional to t*) is the optical depth along the observer’s line of sight at direction cosine /j and 
radial coordinate r. The constant of proportionality (which is, itself, proportional to the 
emission measure) will not be determined in this work. Equation (4) must be solved numer¬ 
ically, except in the case of (3 — 0 (constant velocity wind). And even so, we only obtain 
solutions for integer values of f3. 

This model assumes implicitly that the sites of X-ray emission are so numerous and well- 
mixed with the primary cool wind component that we can treat the wind as a two-component 
fluid. It also neglects non-radial velocity components, including small scale fluctuations 
like turbulence. Note that the assumption of purely radial velocities is what allows wind 
absorption to break the symmetry of the line, since redshifted emission always arises at higher 
line-of-sight distances (and therefore higher optical depths) than the blueshifted emission (see 
the contour plots in Figures 1 and 2). 

Similar treatments of radiation transport have previously been fitted to the global form 
of low-resolution spectra (Hillier et al. 1993; Feldmeier et al. 1997). MacFarlane et ah (1991) 
and Waldron & Cassinelli (2001) have modeled line profiles from discrete spherical shocks. 
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Elsewhere we report on the investigation of non-spherical models and their applicability to 
hot stars (Kramer et ah 2003; Tonnesen et al. 2003). 

For this study, we have adopted the Owocki & Cohen (2001) model and performed 
fits on eight strong lines in the Chandra HETGS spectrum of the 04 supergiant ( Puppis, 
extracting best-fit values and associated confidence limits for three model parameters: q, R 0 , 
and r*. 


3. Fitting the Model to Observed Line Profiles 

Our data set consists of the ±1 order MEG spectrum from a 67 ks observation of the 04f 
star ( Puppis first reported on by Cassinclli et al. (2001). The FWHM of the MEG spectral 
response is AAmeg = 0.023 A (Chandra X-Ray Center 2001) 4 . All the distinguishable lines 
in this spectrum are many times more broad, allowing their profiles to be well resolved. The 
breadth of the lines means that many of them are contaminated by emission from neighboring 
lines. After eliminating He-like fir triplets as unsuitable for fitting due to excessive blending, 
we identified other potential blends by visual inspection of the spectrum and by referring 
to the line strengths calculated by Mewe, Gronenschild, & van den Oord (1985, Table IV) 
and those in the Astrophysical Plasma Emission Database (APED, Smith et al. 2001) 5 . We 
consider a line with rest wavelength Ao to extend over a wavelength range defined by 

A° ^1 + + AAmeg > A > Ao ^1-— AAmeg • (5) 

The widths of neighboring lines are calculated the same way, and any overlap in the ranges 
is excluded from the fit (see Table 1 for the wavelength range over which each fit was per¬ 
formed). We adopted the terminal velocity value determined by Prinja, Barlow, & Howarth 
(1990), Woe = 2485 km s” 1 . 

We numerically integrate equation (4) over the desired wavelength range using software 
written in Mathematica. The resulting profile is convolved with a Gaussian representing 
instrumental response, binned identically to the data, and normalized to predict the same 
total number of counts as were observed over the same wavelength range. 

To properly treat the statistics of Poisson-distributed, low-photon-count data, we use 
Cash’s C (Cash 1979) as the fit statistic. Fits are performed by calculating C on a grid in 
parameter space. The coordinates of the parameter-space point producing the minimum C 

4 The Chandra Proposers’ Observatory Guide is available at http://cxc.harvard.edu/udocs/docs/ 

5 The Interactive GUIDE for ATOMDB is available at http://obsvis.harvard.edu/WebGUIDE/ 
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value C' m j n are taken as our “best-fit” parameters. We then use limits on AC = C — C m i n to 
define our confidence regions, as described by Cash (1979). The grid is expanded as needed 
until the entire confidence region is encompassed within it. In all our fits we held (3 constant 
at [3—1 (which is very close to the typical O-star value of 0.8, Groenewegen et al. 1989) 
and varied R 0 , q, and t*. 

To confirm the confidence limits derived using the AC statistic, we carried out the 
fitting procedure on monte carlo simulated data sets and compared the parameter-space 
distribution of simulated-data fit parameters to the calculated AC confidence region for 
each line. The monte carlo simulations gave results that were consistent with those given by 
the AC statistic. 

The quality of the fits was evaluated using the Kuiper statistic, a variant of the Kolmogorov- 
Smirnov test. Significance levels were determined using monte carlo techniques. We obtained 
two slightly different distributions of the Kuiper statistic depending on whether the simu¬ 
lated data sets were compared directly to the parent model (giving significance level ao) or 
compared to a new best-fit model found by performing the fitting procedure on the simulated 
data set (giving significance level a RF )- 

In all, we fit eight lines between 6.18 A (Si XIV) and 24.78 A (N VII). The results are 
listed in Table 1, and shown in Figure 1 for two representative lines. All the fits but one are 
formally good according to both distributions. The value of a 0 = 0.039 for the 15.262 A fit 
does not meet the criterion a > 0.05 for a formally-good fit. This line may be contaminated 
by the Fe XIX lines at 15.198 A and 15.3654 A (APED). The relatively-low significance 
values for the N VII line may simply be the result of random variation, or could be a sign 
that there are resolved spectral features not explained by this simple model. In any case, 
the fit is formally good. The 16.787 A fit may also be affected by a blend with an Fe XIX 
line, this one at 16.718 A. This would add flux to the blue edge of the line, increasing its 
skew and explaining the relatively high values of r* and q. 

To demonstrate the typical range of models that can be fit to one line, we show the 
best-fit and two extreme models superimposed on the Fe XVII line at 17.05 A in Figure 2. 
The two extreme models are for the parameter sets that have the largest and smallest values 
of r* within the 95.4% confidence region. 

Finally, we summarize the best-fit and 95.4% confidence limits of the three free model 
parameters for seven of the eight lines in Figure 3 (the values for the Si XIV line are not shown 
because it is at a much shorter wavelength and its fit parameters are poorly constrained). 
Trends in these fitted parameters and their implications are discussed in the next section. 



4. Discussion 


The primary result of the analysis presented in this paper is that the X-ray emission 
lines in the prototypical O supergiant, ( Puppis, can, for the most part, be adequately 
fit with a spherically symmetric wind model having a small number of free parameters. 
Furthermore, the derived parameters are quite reasonable in the context of most wind-shock 
models, being consistent with hot plasma uniformly distributed throughout the wind above 
a moderate onset radius, X-ray emitting plasma extending out to the wind terminal velocity, 
and the need for the inclusion of some wind attenuation. 

In detail, however, some interesting trends emerge. First of all, the amount of wind 
attenuation is significantly smaller than what one might expect from a spherically symmetric 
smooth wind, given what is known about this star’s mass-loss rate and wind opacity. There 
have been various calculations of the wind optical depth (often expressed as the radius of 
optical depth unity) as a function of wavelength for this star. They range from values 
much bigger than what we derive here (7 < r* < 30 calculated by Hillier et al. 1993, using 
M = 5.0 x 10 -6 M 0 yr -1 , i?* = 19 R e , Uqo = 2200 km s' 1 ), to values modestly bigger 
(4 < r* < 8 calculated by Cassinclli et al. 2001, using values from Lamers & Leitherer 1993: 
M = 2.4 x 10' 6 M 0 yr' 1 , A* = 16 R & , Voo = 2200 km s' 1 ). Note that Hillier et al. (1993) find 
different values for Ri depending on whether helium recombines or remains ionized in the 
outer wind, but at energies above 0.5 keV (where all of the lines presented here occur) there 
is little difference between the two scenarios. More recent stellar parameters determined by 
(Puls et al. 1996) (M = 5.9 x 10“ 6 M 0 yr' 1 , f?* = 19 f? 0 , and = 2250 km s' 1 ) agree well 
with the values used by Hillier et al. (1993), but would increase the Cassinelli et al. (2001) 
r* values by a factor of 2, given the same opacity (see eq. 2). 

If we accept the r* values derived from our fits, then the disparity between those values 
and the ones mentioned above suggest that either the mass-loss rates or wind opacities are 
being overestimated in previous calculations. The mass-loss rate of ( Puppis is by now quite 
well established using UV absorption lines and Ha, although improper ionization corrections 
or clumping could lead to systematic errors. The wind opacity determination seems much 
more uncertain, both because of the inconsistent values in the literature and because of the 
difficulty in determining the ionization state of the wind (MacFarlane et al. 1993; MacFarlane, 
Cohen, & Wang 1994). Recent advances in stellar atmosphere modeling may help to improve 
these determinations (Pauldrach, Hoffmann, & Lennon 2001). Another means of lowering 
the wind attenuation is to clump the wind into small clouds that are individually optically 
thick rendering the wind porous and enhancing the escape probability of X-ray photons, thus 
lowering the mean wind opacity. This would also affect the mass-loss rate diagnostics, but 
is, itself, an independent effect. 



An even more curious result of the r* fits is that they are nearly independent of wave¬ 
length. This is surprising because photoionization cross sections should scale roughly as a 
power of wavelength between A 2 and A 3 (Hillier et al. 1993). It is possible that the distri¬ 
bution of ionization edges could conspire to make this relationship much flatter over a small 
range of wavelengths (as the calculations from Cassinelli et al. (2001) seem to indicate). But 
wind clumping might play some role, here too. If the wind opacity is dominated by clumps 
that are individually optically thick across the wavelength range, then the opacity ceases to 
be a function of wavelength and instead depends on the physical cross sections of the clumps 
themselves. We note that the UV line opacity necessary to explain the observed absorption 
line profiles could, in principle, still be provided by the tenuous inter-clump wind, as the line 
cross sections are much bigger than the X-ray photoionization cross sections. 

If R\ 3> R 0 , the line profile is insensitive to changes in R 01 since emission much below 
Ri is largely absorbed by the wind. The values we find for R 0 , though generally small, are 
comparable to our values of Ri (from r* by eq. 3). It is hard to assess these relatively small 
onset radii in the context of the small (sometimes surprisingly small) values claimed on the 
basis of observed f ji ratios in He-like ions (Cassinelli et al. 2001; Kahn et al. 2001; Waldron 
& Cassinelli 2001). This is partially because we do not fit the profiles of any He-like lines 
(they are too blended) and partly because the most extreme results (smallest value for Ri) 
are for S XV, which is a higher ionization stage than any of the lines we fit. 

The fit results for the parameter q indicate that there is not a strong radial trend in 
the filling factor. One might expect some competition in a wind shock model between the 
tendency to have more and stronger shocks near the star, where the wind is still accelerating, 
and the tendency for shock heated gas to cool less efficiently in the far wind, where densities 
are low. Perhaps these two effects cancel to give the observed q ~ 0 relationship. 

In conclusion, the simple, spherically symmetric wind shock model is remarkably con¬ 
sistent with the observed line profiles in ( Puppis, providing the most direct evidence yet 
that some type of wind-shock model applies to this hot star. However, there are indications 
that the absorption properties of the wind of ( Puppis, and perhaps other hot stars, must 
be reconsidered. We will fit this same model to the Chandra spectra of other hot stars in 
the future. But the lack of strong line asymmetries in stars such as ( Ori and 5 Ori and the 
narrow lines in 6 1 Ori C and r Sco indicate that spherically symmetric wind-shock models 
with absorption may not fit the data from these stars as well as they do ( Puppis. 
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Tabic 1. Best-Fit Parameters with 95.4% Confidence Limits 


Ion 

(A) 

A 0 /AA 

2-Uoo/An 

q 

R 0 

n 

*^min 
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cio 

«RF 

N VII 
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_n ^+ 0 - 6 
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o +1.3 
z ’-0.6 

0 5 +2 - 
u -°-0.5 

-l.ii 

l.ii 

92 
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0 VIII 
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14 

—0 i +0 - 6 
U ’ i -0.4 

1 2" 1 " 2 ' 1 
i ’ z -0.2 

9 c+2.5 

-1.04 

1.14 

69 

0.941 

0.790 

Fe XVII 

17.054 

741 

12 

—0 6 +0,4 
U.O_0.2 

1 4+°- 6 

0 5 +1 ' 

u - J -0.5 

-0.74 

1.16 

54 

0.971 

0.950 

Fe XVII a 

16.787 

730 

12 

0.4lg;| 

l.ol 2 - 3 

4-5+I;I 

-1.16 

0.74 

53 

0.517 

0.490 

Fe XVII a 

15.262 

664 

11 

-0.818:1 

1 4+ 1 - 1 

i-5l?:i 

-0.81 

1.19 

50 

0.039 

0.151 

Fe XVII 

15.013 

653 

11 

—0 2 +0,4 
u - z -0.3 

i 4+0- 6 

l-Ol 0 .5 

-1.19 

0.79 

49 

0.887 

0.800 

Ne X 

12.13 

527 

9 

—n 4+ 0 - 5 
U ’^-0.3 

1 4+ 0 - 6 

l.Ol}; 5 

-1.23 

0.03 

26 

0.432 

0.100 

Si XIV b 

6.18 

269 

4 

-o.2t 0 ; s 

1-41l4 

i c;+ 5 - 5 

-0.98 

0.59 

16 

0.857 

0.622 


Note. — The width of the instrumental response in wavelength units is AA = AAmeg — 0.023 A, or 
in velocity units Av = cAX/Xq. The scaled wavelength x = (c/v 00 )(X — Xo)/Xq. N q b s is the number of 
wavelength bins included in the fit. 

a The anomalous results for these lines may be due to contamination. 
b At the 95.4% confidence level upper limit, q is unconstrained for this fit. 
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Fig. 1.— Two representative lines with best-fit models. Shown are Chandra MEG spectra 
of (a) Fe XVII at 15.01 A and (c) O VIII at 18.97 A (in gray), with corresponding best-fit 
profiles (in black). Laboratory rest wavelengths are indicated by the vertical dashed lines. 
To the right we show (b and d) contour plots representing the models. The observer is 
located at (p/R* = 0, z/R * = oo). The inner circle is of radius i?*, the outer circle of radius 
R 0 . Gray contours are curves of constant line-of-sight velocity component in units of Voq. 
Dashed contours are curves of constant optical depth (r = 0.5,1,2), integrated along the 
line of sight. 




Fig. 2.— Models at extremes of the confidence region. Shown are the Chandra MEG 
spectrum of the 17.054 A line of Fe XVII (in gray), and (in black) the best-fit model (r* = 
0.5), the fit with r* held at its 95.4% confidence upper limit (r* = 1.5) and the fit with r* 
held at its 95.4% confidence lower limit (r* = 0.0). The contour plots are the same style as 
in Figure 1 and correspond to (b) the r* = 0.0 model, (c) the best-fit model, and (d) the 
r* = 1.5 model. 
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Fig. 3.— The best-fit values of (a) q, (b) R a in units of i?*, and (c) r* for seven of the eight 
lines (bullets), along with the range given by the 95.4% confidence limits (triangles). See 
Table 1. 






